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ABSTRACT 


The  report  analyzes  the  two-phase  thrust  bearing  near  the 
axis  of  rotation  where  the  radial  pressure  gradient  is  non 
zero  and  the  viscosity  is  negligible  for  the  gaseous  phase. 
The  analysis  agrees  with  the  previous  work  by  Sparrow  and 
Gregg  (2)  for  condensate  thickness  for  the  limiting  case  of 
zero  pressure  gradient.  With  radial  pressure  gradients, 
the  condensate  thickness  is  seen  to  have  the  chance  of  re¬ 
evaporation  due  to  the  change  of  vapor  temperature  in  the 
radial  direction.  The  sensitivity  of  the  condensate  thick¬ 
ness  to  temperature  difference  is  noted.  This  work  is  part 
of  the  complete  analysis  of  the  two-phase  thrust  bearing 
being  studied. 
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1.0  INTRODUCTION 


The  study  presented  in  this  report  contains,  the  analysis  of  the  two 
phase,  externally  pressurized,  thrust  bearing  in  the  area  near  the  axis 
of  rotation.  This  area  of  the  thrust  bearing  is  considered  to  be  Regions 
II  and  III,  as  indicated  previously  in  Reference  1.  (See  Figure  1.) 

Sparrow  and  Gregg  (2)  solved  the  problem  of  condensation  on  a  rotating 
disk  with  no  radial  pressure  gradient  by  using  the  similarity  sub¬ 
stitutions  which  Von  Karman  previously  introduced  for  the  rotating  disk 
problem,  that  is. 

This  is  possible  because  the  temperature  boundary  conditions  at  the  inter¬ 
face  and  disk  surface  are  independent  of  the  radius  r.  The  substitutions 
reduce  equations  (1)  and  (2)  to  differential  equations  in  z  only.  The 
condensate  rate  dependent  on  w  is  then  seen  to  satisfy  the  temperature 
boundary  conditions,  as  well  as  the  equations  of  motion. 

Sparrow  and  Gregg  solved  the  set  of  nonlinear  differential  equations  by 
computer,  and  noted  that,  if  in  equations  (1)  and  (2)  after  the  similar¬ 
ity  substitutions  are  made  the  inertial  terms  are  neglected  compared  to 
the  viscous  forces,  the  solution  is  easily  evaluated.  The  results  of  this 
formulation  led  to  the  results 

The  same  problem  can  also  be  analyzed  by  using  the  approximate  integral 
method.  (See  Appendix.)  This  method  again  leads  exactly  to  equation  (a) 

for  the  condition  that  /  C A7~ ]  ^  is  small.  The  range  of  interest  of  this 

parameter  for  the  problem  under  study  can  thus  be  considered  equivalent 
to  neglecting  the  inertial  terms .  The  results  of  this  analysis  are  shown 
in  Figures  2  and  3,  and  explained  in  the  Appendix. 


(a) 
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The  corresponding  problem  of  a  thrust  bearing  does  not  generally  permit 
the  assumptions  of  zero  radial  pressure  gradient.  The  acceleration  of 
the  vapor  flow  in  the  radial  direction  does,  in  fact,  have  a  pressure 
gradient  in  the  case  of  a  bearing,  thus  requiring  the  consideration  of 
the  pressure  gradient  and  its  effect  on  the  condensate  thickness  and 
temperature  in  the  vapor . 


2.0  ANALYSIS 


The  problem  can  be  conveniently  formulated  in  cylindrical  coordinates, 
where  the  coordinate  of  z  is  taken  as  the  axis  of  rotation.  Angular 
symmetry  is  assumed,  that  is,  b/dQ  =  0.  The  physical  problem  is 
represented  by  a  rotating  disk  upon  which  condensate  forms  due  to  a 
temperature  difference  between  the  surface  of  the  disk  and  vapor. 


Cross  Section 
of  Condensate 


6>  kJ 


The  analysis  of  the  two  regions,  that  is,  II  and  III,  is  done 
simultaneously  with  the  following  assumptions: 

1.  Interface  shear  at  h(r)  between  the  vapor  and  liquid 
phases  is  assumed  zero;  hence,  ='^=0atz=h. 


2.  The  radial  pressure  gradient,-^,  is  the  same  for  the 

liquid  and  the  vapor  (that  is,  ^  =  0)  and  for  this 

problem  will  be  an  assumed  known  function  determined 
solely  by  the  vapor  flow  and  geometry  of  the  region. 
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2.1  Formulation  of  the  Equations  of 
Motion,  Continuity,  and  Energy 

The  equations  of  motion  for  an  incompressible  fluid  in  cylindrical 


coordinates  (3),  with  axial  synmetry,  can  be  written  as 


-TC  -  '  '9^ 


^*2^ 4  ^  nA-V"  _ 


'  [ 


t- 


T)A>'‘  VU  A/  /  J 

X  I 

i'V  '  J 


(1) 


(2) 


t 


“ 


(3) 


The  equation  of  continuity  if 


The  energy  equation  can  be  written  as 

(5) 

which  is  the  well  known  Nusselt  approximation  where  dissipation  is 
neglected  and  the  predominate  mode  of  heat  transfer  is  by  conduction 
axially  across  the  condensate  thickness. 


At  this  point,  it  is  well  to  consider  a  limiting  condition  of  the 
equations  of  motion  and  its  effect  on  the  formulation  of  the  method 
of  solution. 


This  problem  does  not  permit  the  similarity  substitutions.  However, 

it  can  be  formulated  by  using  the  momentun  integral  technique.  The 

assumptions  that  "SliT >'2^itr  are  made  in  addition 

to  the  previous  assumptions.  Note  that  as  ^2 -»»  ,  due  to  symmetry,  2^  a 

/  ^ 

if  no  sources  are  present.  Furthermore,  if  the  equations  of 

'D'V 

motion  will  permit  the  similar  solutions  of  Von  Karman  (that  is,  etc)  . 

The  special  case  of  A  ■  0,  is  that  which  was  solved  by  Sparrow  and  Gregg. 
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The  problem  of  condensation  with  0  will  be  formulated 

for  the  general  case  where  the  Inertial  terms  are  not  neglected. 
The  value  of  is  assumed  a  known  function  of  r.  However, 

since  it  has  been  shown  in  the  case  for  '2?  O  that  neglecting 
inertial  terms  in  comparison  to  viscous  terms  was  equivalent  to 
the  case  of  small  film  thickness  or  low  heat  transfer  rates, 
this  assumption  will  be  made  subsequently  to  simplify 
the  problem  where  the  pressure  gradient  is  non  zero. 


2.1.1  Momentum  Integral  Equations 
for  the  General  Case _ 

Equations  (1)  and  (2)  ,  neglecting  higher  order 


viscous  terms  (that  is,  ,  ^  ^~^bCr'  ,etc.) 


become 


(6) 


'3^ 


(7) 


From  equation  (4),  w  at  z  =  h,  can  be  written  as 


Equation  .(6^  can  be  integrated  over  the  condensate 
layer  thickness,  resulting  in 


(8) 


(9) 


where  s.  —J. 

V- 


(I'ffV 

til) 

(12) 
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The  following  identities, 

-  X  ®  ^  Cau) 

and 


(13) 


(14) 


M'S. 

are  substituted  into  equation  (9  )  resulting  in 

4, 

ru  Jo  ^  7" 

Similarly,  e^^uation  7,  can  be  integrated,  which  results  in 


(15) 


=■  ?!< 
r 


where 


Zm  -  >/'2ur/ 


The  following  identities  can  be  made: 

^  1/^  *  J  ^ 

o 


(16) 


(17) 


(18) 


(19) 


(20) 


c 
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f  _  Id.  - 

where  7  =  ^  and  ?  =  (22a) ,  (22b) 

Using  equations  (18)  --  (22b),  equations  (15)  and  (16) 
can  be  written  as 


^(9.  ifh  X  (e.  V^) .  Vd^  ( A*-4 j.  .x‘) 

- 

^  f 

and 

Equations  (23)  and  (24)  are  thus  the  Momentum  Integral 
Equations  for  the  case  of  rotational  symmetric  flow. 


(23) 


(24) 


In  order  to  solve  these  equations,  velocity  profiles  for 
both  u  and  v,  (hence,  f  and  g)  are  assumed  to  be  polynomials 
in  z . 


2. 1.1.1  Velocity  Profiles 
Radial  Velocity  Profile 

The  radial  velocity  profile  f.,  is  assumed  to  be  a  fourth 
order  polynomial  in  q,  where  q  -  z/h. 


Hence , 


(25) 


-6- 


The  conditions  imposed  on  the  radial  velocity  profile  are 
JjrOj  ”  (FU)  -t  ^iO^) 

j.K,  u.ro ^  “n  -  -  r ] 


In  terms  of  f,  and  t) 


These  boundary  conditions  are  representative  of  the  physical  con¬ 
ditions  imposed,  as  well  as  satisfying  the  original  differential  equations. 


Using  the  boundary  conditions,  the  constants  a,  b,  c,  d, 
are  evaluated  and  are 


> 

(26a) 

m 

b  =  -  4,,/^ 

C  =  "  Z,  ■*^^2.(^1'!'^) 

(26b) 

(26c) 

d  -  /-  <4,/^  -  54/3 

(26d) 

where  ^  0 

(27) 

-  j-  ( 

(28) 

Finally,  f  becomes 

/-  0-yJ  O’j/-  ^  ^ ) 


(29) 


Tangential  Velocity  Profile 

Similarly,  the  tangential  velocity  profile  can  be  written  as 

^  a.  ■i’ hy-f  *  ep"*  (30) 

The  boundary  conditions  are; 

•j..  ,  K-.rw  ^ 

>  lyf  ''  0<f  A. 


zrdv 

7^ 


or  in  terms  of  g  and  r| 

T’‘, 

r‘‘  r'  hrr(^*^’^) 

The  constants  evaluated  are 


a  =  /— S 

b  = 

(31a) 

(31b) 

c  =  0 

(31c) 

d  =  -2S  - 

(31d) 

®  ■  S4  4/3 

where  ^  C/V"  j 

0\  IV  TlAf) 

(31e) 

Finally,  g  becomes 

'*/-S  -k  y/( yj)Cf4  2t^) 

(32) 
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Following  the  line  of  thinking  similar  to  Bohlen-Holstein 's 
method,  shape  parameters  can  be  introduced  as  in  the 
Karmen- Pohlausen  technique.  (3) 


It  will  be  apparent  later  that  it  was  convenient  to  intro¬ 
duce  the  so-called  "shape  parameters",  which  in  this  case 


are  defined  as 

"T  ^ 

(33) 

Ki-- 

(34) 

K,  -  A*  jr 

^  A 

(35) 

0 

(36) 

Also,  let 

(37) 

Hence,  K,  *  V 

(38) 

zr 

II 

(39) 

A 

A. 

(40) 

(41) 
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From  the  definitions  of 


and  S,  which  are 


<(>1^  hLEOnK  t^roJ- 


•i  XT 

•i  "XJ 

(42) 

<k  .-  F('0  ^ 

O  IT 

ij  vl'  _  k}-(in 

i)  MJ-  i> 

(43) 

'k^ 

id:^dv 

'  i  A, 

<>  V  (tA, 

(44) 

and  S  = 

(45) 

aQd  defining  Xj  *ik  j  “  ik^-ST 


and  ^4.  -  ^2^  (1^ j and  S  can  be  written  in 


terms  of  ^  ,  T-j 

and  'X  ^  alone . 

hence : 

4^1  •  2:1 

Ai 

(46) 

4^2  * 

^  M  d  ^4-  -  d 

T4  ^  Jli 

(47) 

4}  - 

f^Xv)  —  ^  ^  g/  ^ 

^A4  CtAy 

(48) 

/-  2jr 

Ai- 

(49) 

This  is  evident  by  the  following  relations: 


-40- 


(50) 


and  similarly, 


From  the  definitions,  !%  it  is  readily  seen 

that  after  integration 


MS  ‘ffO  ^072 


t  si.i, 

iUtv  "^fTTl 


.  I 


ysli 


-  2  3^zS  ^ 

ao2^  ?o  aozf 


(51) 


(52) 

(53) 

(54) 


-  5*  4,  19 

iBi»H  Xisro 


(55) 


Hence,  it  is  seen  that  the  shape  factors  satisfy  the 


universal  relationships,  that  is 

K.  •  FJ^->  , 


9  XT 


?  XT 


(56) 


(57) 


(58) 


(59) 
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Having  defined  the  shape  factors  their  expression  can  be 
introduced  into  the  momentum  integral  equations. 


(60) 


Equation  (23)  can  be  written  as 

+  r 

Multiplying  equation  (60)  by 
and  recalling  the  definrtions  of 

given  by  equations  (52),  (53),  (54),  and  (55),  and  the  relation 

=  F,CK,) 

yir  ^  ^ 

one  obtains,  using  the  definition  of 
^  hence,  the  following  equation 

2£iS  ^  -!<^h  I 

i  ^  Ks  ^  Kl  Tl  3 

V 


(61) 


(62) 


Let 


(63) 


and  F^(K(;')  = 

ii. 


(64) 


one  obtains 


(65*) 


-12- 


In  a  similar  manner,  the  second  momentum  integral  equation  can 


be  written  as 


where 


^  .  fT 

2.  ^  ' 


'£  •  f<\(  tL  ^  ) 

^  (  tl.  -.  L  ) 


from  the  definition  *^2’^“ 


one  can  obtain  d.  A  ^2.)^  'y'd}r-nrcii^ 

Equation  (65)  and  (66)  subject  to  the  conditions  of  equation  (70) 
or  (71)  contain  three  unknowns;  that  is,  U,  V,  and  h. 

Equation  (71)  is  merely  an  identity  introduced  for  convenience  in 
performing  the  numerical  solution.  The  third  equation  needed  for 
the  complete  solution  is  one  which  is  obtained  by  considering  the 
energy  equation  which  relates  the  value  h  to  the  value  U  under  the 
temperature  differentials  imposed. 


(66) 

(67) 

(68) 

(69) 

(70) 

(71) 


The  final  set  of  differential  equations  will  appear  as 


J-cr 

=  V  (  0;  V .  ■  ) 

(72a) 

iX  = 

'V  ( 1^, 

dn. 

(72b) 

=■ 

dAj 

k  [0,  Ar  -) 

(72c) 

where  (72c)  Is  obtained  from  energy  considerations. 
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These  equations  generalized  to 'include  energy  equations  correspond 
Karmen-Pohlausen  method,  but  in  rotational  symmetric  coordinates. 


2 . 2  Special  Case:  Viscous  Forces  Predominant 

As  mentioned  before,  the  range  of  interest  in  the  problem  would 
be  that  corresponding  to  the  equations  with  no  inertial  terms. 
If  this  is  assumed,  the  following  relationships  are  known: 

V  =  ru) 

S  =  0  j 

,  ,  -  ^  XT  'LT  ' 

g  -  I  ^  --TTf- 

From  equation  (61) 


Tin, 

T 


Z^V 

1^ 


But  also  from  equation  (60),  neglecting  the  inertial  terms. 


7^,74, 

Since  4s  2,  Che  function  £  becomes 


to  the 


(73) 

(74) 

(75) 

(76) 

(77) 

(78) 
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Energy  Equation 

Using  the  Nusselt  approximation,  the  energy  equation  is  reduced 
to  the  Fourier  conduction  law  across  the  condensate  thickness . 
Hence,  in  the  annular  region  of  area  2irrdr,  the  heat  flow  dq  is 


(79) 


This  exchange  of  heat  is  related  to  the  amount  of  condensate 
formed  by  the  relationship 

(i%~  ^  (3,  (80) 

where  is  the  amount  of  condensate  added  in  the  differential 
area  under  consideration. 


The  flow  of  condensate  in  the  radial  direction  is  given  by 

(81) 

whereiTis  the  average  velocity,  radially. 


Using  equations  (78)  and  (75),  one  obtains 

(5  s  nejL  ^ )[  ^ 

Defining  ^CA')  »  AvJ,  I-I  s  ^  ^ 

AW**-  ^ 

(5  jr  V t  4  ) 

Therefore,  since  Q  =  9-  (H> 

do 


(82) 

(83) 

(84) 

(85) 

(86) 
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Using  equations  (84),  (P5.*  {86)  and  {79) 

-h  4  2/l^/-^^)J  (87) 

3 

Making  the  following  substitutions 

=*  m  r-7W  45*  j-  C 

^  x^/w ' 

reduces  equation  (87)  into  the  following  form. 


The.  Clausius-Clapeyton  equation  relates  the  pressure  to  the  temperature, 

^  .  2MP 

^r  (89) 

Integrating , 


(90) 


Hence , 


f(^) 


1+  lyis  jtwu£. 
AT  AM  Br 

I- 

XM  15 


Finally,  one  can  let  Y-(«V>/4©‘) 


(91) 

(92) 


and  obtain 


4^=-6  J1  -  i  i  ? 

^  3a  i  I  rrXTT  i 


(93) 


All  that  is  needed  is  to  prescribe  V(a) ,  and  solve  equation  (93) 
numerically.  Note  chat  at  R  -  0^  — !►  ^  V* ■ 
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Particular  Case;  (P)r)  Is  Assumed  Parabolic 
A  parabolic  pressure  distribution  Is  assumed  where 
The  boundary  conditions  are: 


^  A*  O 


P 


a  +  br 


2 


Hence , 


!  +(^.,)  |2. 


where 


oi  s.  Px 

?S 


From  equation  (96) , 


At 


(94) 

(95a) 

(95b) 

(96) 

(97) 


But,  iU)  - 

r 

Hence,  s  2  -61  «  CqiOW. 

P  /lo'y*-  ’ 


Using  equation  (96),  (97)  and  (98)  results  in 

^  =L  I  ill)  J-  -e  'iji) 

^  ^  S  IH  5  s 

tL%  ')m  l-h  ^  ”0  ^ 

I- 


(98) 


(99) 


(100) 
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Equation  (99)  can  be  simplified  further  by  letting^T/V ^^/y/^resulting  in 

-fjr  (101) 

subject  to  the  initial  condition  U(0)  =  1. 

Equation  (101)  can  be  solved  numerically  by  using  a  Runga-Kutta. 
integration  procedure. 


APPROXIMATE  SOLUTION 

Equation  (101)  possesses  a  solution  which  by  general  theory  of  first  order 
equations  can  be  written  as 


t/TE)e  =  J  ^ 


where  p/ »  6 

3J 


(102) 

(103) 

(104) 


r(i) 

The  integration  of  is  not  possible  analytically,  due  to  (iCV. 

However,  under  certain  conditions ^(^^may  be  simplified  in  the  following 
manner : 


The  value  of  the  function  In  (1-  t^)  can  be  written  as 


for 


JUO-y) 

/  '  _«v  -  «' 


(105) 


yj  -at  — 

as  compared  to  the  exact  value  of  0.2225. 
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If  ^  (l“^J ,  since  the  maximum  value  of  ^2/  ,  the  approximate 
polynomial  expression  will  be  good  for  the  full  range  of^for 
For  values  of  ,  the  accuracy, "of  course,  decreases.  However, 

for  estimation  purposes,  it  is  desirable  to  formulate  the  solution  to 
equation  (101)  analytically. 


Using  the  approximation  of  equation  (105),  but  omitting  the  cubic  term, 
transform  equation  (101)  to 


where  Ks.  *  Ik  .  K,  = 


(106) 

(107) 


Since  <i  < I  ,  the  denomtuator  of  equation  (106)  can  be 
written  as  t ^ ^  where  ^  / 

Hence,  [[' 1 []  I' l<'|  ^  0  ^ 


(108) 


neglecting  higher  order  terms  , 

Hence,  I- ^  (109) 

Using  the  results  of  equation  (109)  in  equation  (102) 

(J  •  I / (110) 

Hence  XT  after  integration  becomes 


>», , 

Figure  /7  shows  a  plot  of  with  ^  as  the  abscessa  and  0(  as  a 

parameter . 


(Ill) 
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3.0  DISCUSSION 


The  dimensionless  condensate  thickness  for  the  special  case  of  zero  pressure 
gradient  is  shown  in  Figure  (3).  The  results  agree  extremely  well  with  those 
of  Sparrow  and  Gregg  (2).  The  distributions  of  condensate  thickness  for 
various  pressure  gradients  resulting  from  an  assumed  parabolic  pressure 
distribution  are  shown  in  Figures  4  to  15. 


The  results  are  divided  into  four  sets  at  a  constant  pressure  ratio  and 
varying  source  pressure.  The  abscissa  of  the  curves  is  the  dimensionless 
radius  ratio  ^  .  The  ordinate  of  the  curves  is  essentially  the  ratio  of 
dimensionless  condensate  thickness  with  pressure  gradient  to  that  of  zero 
pressure  gradient,  that  is, 


T  Jj_  /  / ZFsC'-o^)  \ 

where  ^  -  A  ^  ^ 


(112) 


i  I 

The  curves  of  (J  versus  the  radius  ratio- ^  ,  show  for  a  given  value  of 
source  pressure  and  pressure  ratio,  L/*^'^increases  with  increasing  tempera¬ 
ture  difference  and  decreasing  radius.  As  the  radius  ratio  approaches  zero, 
the  value  of  U  ^approaches  1,  independent  of  the  temperature  difference  AT"  • 

As  the  value  of  the  angular  velocity <0  ,  approaches  zero,  the  value  of  U  ^ 
becomes 


0 


'4 


a  K. 


3^)4  J 


^4. 
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where  K  is  determined  as  a  function  of  J*,  and  AT.  For  values  of  large « 
I  .  Therefore, 


3  ^ 

(i-<)  i 


'4 


An  interesting  aspect  of  the  condensation  in  the  central  region  of  the 
bearing  can  be  seen  by  referring  to  Figure  18.  Figure  18  shows  a  plot 
of  the  dimensionless  condensate  thickness  iJ  ^versus  the  temperature  dif¬ 
ference  7^  -TJv  at  5  =  1  for  various  pressure  ratios  and  at  a  constant 
source  pressure  of  200  psi.  From  this  curve,  it  is  readily  seen  that  the 
height  of  the  condensate  layer  at  the  outer  radius  of  the  central  region 
is  very  sensitive  to  the  temperature  difference.  For  a  given  pressure 
ratio  o<'  ,  a  narrow  range  of  47~exists  at  the  extremes  of  which  the  bearing 
will  respectively  be  dry  or  attain  a  maximum  thickness  of  condensate  at 
the  exit  of  this  region.  This  effect  of  temperature  difference  on  the  con¬ 
densate  layer  can  thus  cause  extreme  variation  in  the  operating  character¬ 
istics  of  a  particular  bearing. 


Figures  4  through  19,  also  indicate  that  for  a  given  source  pressure  and 
pressure  ratio  there  exists  a  temperature  difference  such  that  the  conden¬ 
sate  layer  becomes  zero  at  ?  =1.  These  points  are  plotted  in  Figure  16, 

7v  -- 

where  the  ordinate  is  chosen  at  ii'stead  the  abscissa  being  oL  , 

the  pressure  ratio.  At  a  constant  source  pressure,  hence  a  constant  source 

temperature,  as  and  hence  a  dry  condition  is  prevalent 

i-  ' 

at  JS  =  !•  As  Ts -Tyn  gets  larger,  a  fully  wet  bearing  at  f  =  1  can  result. 
Hence,  for  any  constant  source  pressure  ,  Figure  16,  shows  that  if 
is  above  the  curve,  a  dry  bearing  will  always  exist  at  §  =1,  the  reverse 
being  true,  if  the  ordinate  is  below  the  curve,  • 
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The  effect  of  L)  versus  ^  in  Figures  4  through  15  shows  that  for  a 
given  4 r >  especially  where  this  ^T'leads  to  a  dry  bearing  at  ^  less  than 
1,  that  the  values  of  0^^^  change'  very  rapidly  near  the  point  where 
reevaporation  occurs.  An  explanation  of  this  fairly  rapid  change  in  con¬ 
densate  thickness  can  be  attributed  to  the  effect  of  condensate  "thin¬ 
ning"  caused  by  a  temperature  difference  which  is  now  negative  (that  is, 
the  disk  temperature  is  higher  than  the  vapor  at  the  particular  radius) 
and  the  continuously  applied  centrifugal  force.  Before  the  temperature 
reversal,  the,se  two  effects  opposed  each  other,  that  is,  thinning  down 
of  the  condensate  increases  the  condensation  rate.  When  this  point  is 
reached,  these  two  effects  are  additive  in  the  sense  that  they  both  work 
in  the  same  direction. 

The  condensation  process  as  formulated  in  the  report  proceeds  until  the 
point  is  reached  where  the  saturation  temperature  of  the  vapor  is 
equal  to  the  given  wall  temperature.  Beyond  this  point  the  validity 
of  the  present  analysis  is  in  question.  The  present  equations  predict 
reevaporation  over  a  small  radial  increment  which  also  infers  a  ne¬ 
gative  temperature  difference  between  the  wall  and  interface.  This 
situation  that  occurs  where  the  liquid  temperature  at  the  wall  is 
greater  than  the  vapor  temperature  is  in  disagreement  with  the  initial 
assumptions  of  thermodynamic  equillibrium.  A  new  self  consistent 
physical  model  needs  to  be  formulated  for  the  evaporation  process. 

Whereas  the  present  analysis  cannot  be  depended  on  to  describe  the 
evaporation  process  precisely,  it  does  suggest  that  the  evaporation 
takes  place  in  a  small  radial  distance.  The  existence  of  a  superheated 
fluid  at  a  temperature  above  the  corresponding  saturation  temperature 
is  a  phenomenon  which  is  attributed  to  a  process  which  is  not  in  thermo¬ 
dynamic  equillibrium  but  rather  one  where  a  boiling  process  is  prevalent. 
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4.0  CONCLUSIONS 


1.  For  a  given  degree  of  subcooling,  the  condensate  film  thickness 
near  the  axis  of  rotation  is  strongly  dependent  on  the  magnitude  of  the 
radial  pressure  drop  (See  Fig.  18).  Hence,  the  ratio  of  mass  flow  be¬ 
tween  the  liquid  and  the  vapor  phases  entering  the  bearing  region  is  in¬ 
fluenced  by  both  these  parameters. 

2.  The  limiting  case  of  zero  radial  pressure  gradient  is  obtained 
from  the  present  analysis  which  agrees  with  the  work  of  Sparrow  and 
Gregg  (2). 

3.  The  limiting  case  of  zero  angular  velocity  reduces  to  the  case 
of  condensation  in  axisymmetric  flow.  The  results  for  this  case  can 
readily  be  obtained  by  using  the  same  curve  in  this  report  as  pointed 
out  in  the  Discussion. 

5.0  RECOMMENDATIONS 

1.  The  results  contained  in  this  analysis  can  be  used  as  the  initial 
conditions  for  the  two  phase  lubrication  equations  which  are  derived  in 
reference  1. 

2.  The  present  analysis  is  applicable  to  the  rotating  condenser. 
Further  study  of  the  Implications  of  the  analysis  to  a  condenser  should  be 
evaluated. 
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NOMENCLATURE 


(All  dimensional  quantities  are  In  consistent  units) 


C  Specific  Heat 

h  Condensate  Thickness 

h*  Displacement  Thickness  of  Condensate,  Defined  by  Equation  (52) 
k  Thermal  Conductivity 

q  Heat  Flow 

r  Radial  Coordinate 

u  Radial  Velocity 

V  Tangential  Velocity 

w  Axial  Velocity 

z  Axial  Coordinate  i 

H  Dimensionless  condensate  thickness,  ^ 

M  Molecular  Weight 

P  Pressure 

P  Source  Pressure 

s 

P  Final  Pressure 

o 

P  Prandtl  Number 

r 

Q  Mass  Flow 

R  Gas  Coolant 

T  Temperature 

T^  Source  Temperature 

T^  Final  Temperature 

T^  Disk  Temperature 

&  Dimensionless  Temperature 
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p 

o 

Pressure  Ratio,  — 
s 

Dimensionless  Thickness, 
Dimensionless  Parameter, 


c  (rs-Tv)/)^P, 


Density 

Dimensionless  Radius, 


Kinematic  Viscosity 


a)  Angular  Velocity 
"X  Heat  of  Vaporization 
^  Shearing  Stress 
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APPENDIX 


The  case  of  condensation  on  a  rotating  isothermal  surface  can  be  analyzed 
using  the  approximate  integral  method,  which  leads  to  a  solution  which  is 
in  agreement  with  that  by  Sparrow  and  Gregg.  (2) 


Since  the  pressure  gradient  is  zero,  (that  is,  F(r)  =  0),  equation  (15) 
can  be  written  as  . 


(T-1) 

and  equation  (16)  as 

(1-2) 

Defining  being  the  condensate  height,  (in  this  case,  a 

constant) , 

and  letting 

u  =A  r'/tj  •= 

(1-3) 

V  = 

(1-4) 

U  = 

(1-5) 

V  =  ol>w<0 

(1-6) 

Equation  (I-l)  becomes 

/  3  Vx  -  v/,  -i 

'  O  * 

and  equation  (1-2) 


(1-7) 


(1-8) 
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The  velocity  profiles  F  i^)  and  G  t'S)  are  assumed  to  be  polynomials,  which 
satisfy  the  boundary  conditions  and  the  original  differential  equations. 

Assume  ,  to  be  a  fourth  order  polynomial,  that  is, 

r(r)=  6^-h6t 

The  conditions’. which  govern  the  constants  a,  b,  c,  d  and  e  are: 

F  (0)  =0 

F'(0)  =  0 

F'(l)  = 

F"(l)  =0 

F  "  '  (0)  =  - 

Thus, 

S  ^  fgi)  ^ 

The  tangential  velocity  profile  G(S)  is  determined  by 
G^('S)  =  a  +  b5+  cjVi/s’  subject  to  the  conditions 
G(0)  s  U» 

G(l)  =  CLUJ 
G'(l)  =  0 

G"(0)  =  0 

ol^S)  =  (1.13) 

Substituting  equations  (I-ll)  and  (1-13),  and  their  corresponding 
derivatives  into  equations  (1-7)  and  (1-8)  and  performing  the  corresponding 
intergrations ,  result  in 


(1-9) 

(I-lOa) 

(1-lOb) 

(I-lOc) 

(MOd) 

(I-lOe) 

(I-ll) 

(I- 12a) 
(I-12b) 
(I-12c) 
(I-12d) 
(I-12e) 


where  //  *  O  iy 
0 


(1-14) 

(1-15) 
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(1-17) 


The  total  mass  flow  in  the  radial  direction  is  given  by 
but  /Lp'Ct) 


where 


Substituting  equation  fl-17)  into  equation  (1-16) 

Since  <i6f  ~ 
and  ^  TC  ifC^) 


(1-18) 

(1-19) 


(1-20) 


Evaluating  equation  (1-20)  in  terms  of  H, 

where  ^  ^  ^T**  (1-22) 

A75L 

The  aolution  of  H,  «(,  and  ^  as  a  function  of  ^  are  sought  which  can  be 
determined  uniquely  from  equations  (1-14),  (1-15),  and  (1-21). 


Simplified  Solution 

Equations  (1-14)  and  (1-15)  are  considerably  simplified  if  H  is 
considered  small,  compared  to  unity. 

Equation  (1-14)  reduces  to 

5S  »5 

and  equation  (1-15)  reduces  to 


(1-23) 


(1-24) 
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This  results  in  ^  ^  ui-/ 

/6  »  ^  /-L^ 

(?-y* 

Substituting  these  results  into  equation  (1-21) 

H .  ij/J'*-  e 

/»  -  CVeJ  © 

•*- 

The  results  agree  with  the  work  of  Sparrow  and  Gregg,  as  shown  by 
equation  (6)  of  this  report. 

A  solution  for  large  values  of  H  was  performed  numerically  by 
solving  for and ^  from  equations  (1-14)  and  (1-15)  for  values 
of  H.  The  value  of  |^(H)  from  equations  (1-14)  and  (1-15)  were 
solved  graphically  with  equation  (I-^  which  involved  versus  H 
for  various  values  of .  This  solution  is  shown  in  Figure  2. 

The  intersection  points  were  cross-plotted  to  obtain  Figure  3, 
whereoC ,  ^  ,  and  H  are  shown  as  functions.of  ^  . 


(1-25), 

(1-26) 


(1-27) 

(1-28) 

(1-29) 


-30- 


FIGURES 


FIO.  2  NUMERICAL  SOLUTION  OF  EOUATIONSa-N.I-IS)  FOR  m,fi  VERSUS  H  AND  EQUATION  (Z-21) 
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FIG.  4  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 
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FIG.  5  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 
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FIG.  6  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 
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FIG.  7  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 
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FIG.  8  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 


FIG.  9  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 


frA‘SS3N>IOIHl  3iVSN3QNOO  SS31N0ISN3Wia 
T 


DIMENSIONLESS  RADIUS  RATIO,  { 

FIG.IO  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 
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FIG.  II  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 
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FIG.  12  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 


FIG.  13  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 


FIG.  14  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 
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FIG.I5  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  DIMENSIONLESS  RADIUS  RATIO 
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FIG.  18  DIMENSIONLESS  CONDENSATE  THICKNESS  vs  TEMPERATURE 

DIFFERENCE  FOR  VARIOUS  PRESSURE  RATIOS  AT  A  CONSTANT 
SOURCE  PRESSURE  OF  200  PSI  AND  f  =  l  O 
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Attn:  Melvin  S.  Day 

Applied  Physics  Laboratory 
Johns  Hopkins  University 
Silver  Spring,  Maryland 
Attn:  George  L.  Selelstad, 

Supr.,  Tech.,  Reports  Group 

Departstent  of  Chemical  Engineering 
Mew  York  University 
Mew  York  53,  New  York 
Attn:  JaiMS  J.  Barker,  Assoc.,  Prof, 
of  Nuclear  Engineering 

Professor  A.  Charnes 
The  Technological  Institute 
Northwestern  University 
Evanston,  Illinois 

Rasldent  Representative 
Office  of  Naval  Reaearch 
c/o  University  of  Pennsylvania 
3438  Walnut  Street 
Philadelphia  4,  Pennsylvania 

Profeasor  P.  R.  Tnimpler 
Towne  School  of  Civil  and  Mechanical 
Baglnecring 

University  of  Pennsylvenla 
Philadelphia,  Pennsylvania 

Jet  Propulsion  Laboratory 
California  Institute  of  Technology 
4800  Oak  Grova  Avenue 
Paeadena,  California 
Attn:  Library 

Illinois  Institute  of  Technology 
Chicago  16,  Illinoia 
Attn:  Professor  L.  M.  Teo 

Profeaaor  N.  C.  Shaw,  Head 
0  Depertuent  of  Nechcnlcel  Engineering 
'  Cemegie  IiuCieute  of  Technology 
pictaburgh  13.  PennaylTenU 


Engineering  Projects  Laboratory 
Haasachuseets  In;iticute  of  Technology 
Cambridge  39,  Masaachuaetta 
Attn:  Dr.  R.  W.  Mann  (Room  3-459  A) 

I 

Haasachusecta  Institute  of  Technology 

Instrumentation  Laboratory 

68  Albany  Street 

Cambridge  39,  Massachusetts 

Attn:  Library,  Wl-109 

1 

Professor  P.  F.  Nartinuiti 
Stevens  Institute  of  Technology 
Hoboken,  New  Jersey 

BaCtelle  Memorial  Institute 
1  505  King  Avenue 

Columbus  1 ,  Ohio 
Attn:  Dr.  Russell  Dayton 

Franklin  Institute 

I  Laboratory  for  Research  and  Development 

Philadelphia,  Pennsylvania 
Attn:  Professor  D.  D.  Fuller 

Library  Inatltute  of  Aerospace  Sclencee 
2  East  64th  Street 
1  New  York,  New  York 

Nr.  G.  B.  Speen 
Sr.  Member,  Technical  Staff 
ITT  Federal  Laboratories 
Division  of  International  Telephone 
1  and  Telegraph  Corporation 

15151  Bledsoe  Street 
San  Fernando,  California 

Aerojet-General  Nucleonics 
P.  0.  Box  86 

I  San  Ramon,  California 

Attn:  Barbara  M.  Frobert 

Aerospace  Corporation 
I  p.  0.  Box  95085 

Loi  Angelea  45,  California 
Attn:  Aeroapeca  Library 

Technical  Reperca  Croup 
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AlResearch  Manufacturing  Company 
Sky  Harbor  Airport 
402  South  36th  Street 
phoenix,  Arizona 
Aten:  Mr.  VanDenwegen 
Librarian 

Mr.  William  D.  Stinml 
Research  Laboratories  Library 
Allls-Chalmers  Manufacturing  Cnpany 
Milwaukee  1,  Wisconsin 

J.  W.  Pecker  Division 
American  Optical  Ccopany 
4709  Baum  Blvd. 

Pittsburgh  13,  Pennsylvania 

American  Society  of  Lubrication 
Engineers 

5  North  Wabash  Avenue 
Chicago  2,  Illinois 

Chairman 

Research  Coosittee  on  Lubrication 
The  American  Society  of  Machanlcal 
Engineers 

United  Engineering  Center 
345  East  47th  Street 
New  York  17,  New  York 

Mr.  James  R.  Kerr,  President 

Lycoming  Division 

AVCO 

Stratford,  Connecticut 

Research  Freda  Ion  Mechanisms 
Division  of  the  Barden  Corporation 
4  Old  Newtown  Road 
Danbury,  Connecticut 
Attn:  B.  L.  Mima,  Vice  Preaidant  - 
General  Manager 

Beemer  Engineering  Cenpany 

Induatrlal  Park 

Fort  Waahington,  Pennsylvania 

Utica  Division 
The  Bandix  Corporation 
211  Seward  Avanua 
Utica,  Maw  York 
Attn:  Mr.  Ruaaell  T.  Dalhith 
guperrlaeiy  toglneer 


Office  of  Assiftant  Director 
(Army  Reactors) 

Division  of  Reactor  Development 
U.  S.  Atomic  Energy  Commission 
Washington  25,  D.  C. 

Attn:  Mr,  Clarence  E.  Miller,  Jr. 

Chief 

Engineering  Development  Branch 
Reactor  Development  Division 
U.  S.  Atomic  Energy  Conmlssion 
Washington  25,  D.  C. 

Headquarters  Library 

U.  S.  Atomic  Energy  Conmlssion 

Washington,  D.  C. 

Chief,  Division  of  Engineering 
Maritime  Administration 
GAO  Building 
Washington  25,  D.  C. 

Cryogenic  Engineering  Laboratory 
National  Bureau  of  Standards 
Boulder,  Colorado 
Attn:  Library 

Mr.  Harold  Hesslng 
National  Aeronautics  and  Space 
Administration 
1512  H.  Street,  N.  W. 

Washington,  D.  C. 

Mr.  Edmund  E.  Bisson 
Chief,  Lubrication  &  Wear  Branch 
Lewis  Research  Center 
National  Aeronautics  and  Space 
Administration 
2100  Brookpark  Road 
Cleveland,  Ohio 

Mr.  Rudolph  E.  Beyer 
George  C.  Marshall  Space  Flight  Center 
National  Aeronautics  and  Space 
Administration 

Guidance  and  Control  Division 
Gyro-Stabilizer  Branch 
Huntsville,  Alabama 

Mr.  H.  W.  Savage 

Oak  Ridge  National  Uboratery 

Poet  Offlca  Box  Y 

Oak  Ridge,  Tannassea 


Bendtx  Aviation  Corporation 
Research  Laboratories  Division 
Southfield,  Michigan 
Attn:  Mr.  Ralph  H.  Larson 

Ml.  C.  R.  Adams 

Physics  Technology  Department 

Aero-Space  Division 

The  Boeing  Company 

Seattle  24,  Washington 

Bryant  Chucking  Grinder  Company 
60  Clinton  Avenue 
Springfield,  Vermont 
Attn:  Mr.  Roald  Cann 

Cadillac  Gage  Cenpany 
P.  0.  Box  3806 
Detroit  5,  Michigan 

Attn:  Mr,  J.  Taylor,  Project  Engineer 

Dr.  Dewey  J.  Sandell 

Director  of  Development 

Carrier  Research  and  Devalopewnt  Co. 

Carrier  Parkway 

Syracuse,  New  York 

Chance  Vought  Corporation 

p,  0.  Box  5907 

Dallas,  Texas 

Attn:  Mr.  R.  C.  Blaylock 

Vice  President,  (Engineering) 

Chryaler  Corporation 
Defense  Operations 
p.  0.  Box  757 
Detroit  31,  Michigan 
Attn:  Hr.  C.  W.  Snider 

J.  D.  Mamarchev  &  Associates 
Consulting  Engineers 
3908  Main  Street 
Houston  2,  Texas 

Mr.  B.  W.  Birmingham 
Cryogenic  Engineering  Laboratory 
National  Bureau  of  Standards 
Boulder,  Colorado 

Curtias  Wright  Corporation 
Wright  Aaronautical  Division 
Department  8332 
Hoed  Ridge,  New  Jersey 

nsjut  bviMTioi 
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Daystroo  Pacific 
9320  Lincoln  Boulevard 
Loe  Angeles  AS,  California 
Attn:  Robert,  H.  Smith 

Special  Project  Engineer 

Ford  Instrument  Company 
31>10  Thomaon  Avenue 
Long  Island  City  1,  New  York 
Attn:  Mr.  Jarvla 

Mr.  Adolf  Egli 
Ford  Motor  Company 
Engineering  and  Research  Staff 
p.  0.  Box  2053 
Dearborn,  Michigan 

Dr.  John  E.  Mayer,  Jr. 

Non-Metalllcs  Section 
Applied  Science  Department 
Scientific  Laboratory 
Ford  Motor  Company 
P.  0.  Box  2053 
Dearborn,  Michigan 

AlResearch  Manufacturing  Division 
The  Garrett  Corporation 
9851  S.  Sepulveda  Boulevard 
Los  Angeles,  California 
Attn:  Jerry  Glaser,  Supervisor 

Mechanical  Lab. -Dept.  93-17 

General  Atomics  Division 
General  Dynamics  Corporation 
P.  0.  Box  608 
San  Diego  12,  California 
Attn:  Mr.  F.  W.  Simpson 

Bearing  and  Lubricant  Center 
General  Engineering  Laboratory 
General  Electric  Company 
1  River  Road 
Schenectady,  New  York 
Attn:  G.  R.  Fox,  Manager 

Mr.  L.  U.  Uinn 

General  Electric  Company 

Aircraft  Accessory  Turbine  Department 

950  Ueetem  Avenue,  Bldg.  3-74 

Lynn,  Meseachusetta 


Research  Laboratories 
General  Motors  Corporation 
General  Motors  Technical  Center 
12  Mile  and  Mound  Roads 
Warren,  Michigan 
Attn:  Mr.  E.  Roland  Makl, 

Mechanical  Development  Dept. 

A.  C.  Spark  Plug  Division 
General  Motors  Corporation 
Milwaukee  1,  Wisconsin 
Attn:  Allen  Knudsen 

Mr.  Walter  Carow 
Kearfott  Division 
General  Precision  Inc, 

1150  McBride  Avenue 
Little  Falls,  New  Jersey 

Grunnan  Aircraft  Engineering  Corp. 
Bethpage,  Long  Island,  New  York 
Attn:  Mr.  David  W,  Craig,  Jr. 

Mechanical  Design  Section 
Engineering  Department 

Hydronautica ,  Incorporated 
200  Monroe  Street 
Rockville,  Maryland 

International  Business  Machine  Corp. 
Research  Laboratory 
San  Jose,  California 
Attn:  Dr,  W,  E.  Unglois 

Mr.  B.  A.  Napier,  Engr.  Div.  Mgr. 
Lear,  Incorporated 
110  Ionia  Avenue  NW 
Grand  Rapids,  Michigan 

Mr.  L.  R.  Barr,  V.  P. 

Lesr-Romec  Division 
Abbe  Road 
Elyria,  Ohio 

Dr.  Calus  G.  Goetsel,  D/S3-30 
Bldg.  201,  Plant  2,  Palo  Alto 
Lockheed  Missiles  &  Space  Co. 

P.  0.  Box  504 
Sunnyvale,  California 


Dr.  J.  S.  Ausiaaa 
Litton  Systems,  Inc. 

5500  Canoga  Avenue 
Woodland  Hills,  California 

Mr.  Don  Moors 
Litton  Systems 
5500  Csnoga  Avsnue 
Woodlsnd  Hills,  Cslifornls 

Mr.  A.  M.  Thomss 
Astro  Division 
The  Marguardt  Corporation 
16555  Saticoy  Street 
Van  Nuya,  California 

Mr.  Kendall  Perkins 
Vice  President  (Engr) 

McDcwnell  Aircreft  Corporation 
Lambert  SC.  -  St.  Louis 
Munieipsl  Airport 
Box  516 

St.  Lolls  3,  Missouri 

Dr.  Beno  Stemlieht 
Mechanical  Technology  Incorporated 
968  Albany-Shaker  Road 
Latham,  New  York 

Mr.  J.  W.  Lower 

Chief,  Engineer-Inertial  Components 
Honeywell  Aero  Divis  ion 
2600  Ridgwsy  Rosd 
Minnespolis,  Minnesota 

Mr.  Carl  P.  Craesser,  Jr. 

Director  of  Reiesrch 
N^  Hampshire  Ball  Bearings,  Inc. 
Peterborough,  New  Hampshire 

Mrs.  Alice  Ward,  Librarian 
Norden  Division  of  United  Aircraft  Corp. 
Helen  Street 
Norwalk,  Connecticut 

Autonetics  Division  (Llbrerlen) 

North  American  Aviation  Inc. 

9150  East  Imperial  Highway 
Downey.  Callfomle 

Northrop  Cotporstion 

Moreir  Division 

1001  Best  Broadway 

Hawthorne,  California 

ACtot  Taehnlcel  Information,  3125 


Mr.  N.  L.  Slomons 
Nortronics 

A  Division  of  Northrop  Corp. 

500  East  Orsngethorpe  Avenue 
Anaheim,  Cslifornls 

Nortronics 

A  Division  of  Northrop  Corporation 
100  Morse  Street 
Norwood,  Massachusetts 
Attn:  Mr.  E.  L.  Swslnson,  Tech.  Aest. 
precision  Products  Department 

Pratt  &  Whitney  Aircreft 
Division  of  UAC  -  CANEL 
P.  0.  Box  611 
Middletown,  Connecticut 
Attn:  Librarian 

Library,  Bldg.  10-2-5 
Radio  Corporation  of  America 
Camden  2,  New  Jersey 

Mr.  Robert  S.  Slegler 
Rocketdyne 

Nucleonics  Subdivision 
6633  Canogs  Avenue 
3  Canoga  Park,  California 

Ryan  Aeronautical  Company 
Attn:  Engineering  Library 
Lindbergh  Field 
San  Diego  12,  California 

Jack  &  Helntt,  A  Division  of 
Tha  Slegler  Corporation 
1725  Eye  Street,  Suite  505 
Washington  6,  D.  C. 

Mr.  Paul  A.  Pitt,  Vice  President 
Engineering  &  Research 
Solar  Aircraft  Company 
2200  Pacific  Highway 
1  San  Diego  12,  California 

Mr.  W.  G.  Wing 

Sperry  Gyroscope  Company 

C-2 

Great  Neck,  New  York 

E.  A.  Babley 

UnloB  Carbida  Muclaar  Company 
Post  Offiea  Box  F 

Oak  Rldga,  Tannaasaa 
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Mrs.  Floranca  Turnbull 
gnglnearlng  Librarian 
Sperry  Gyroscope  Company 
Great  Neck,  New  York 

Sundstrand  Avlatlon-Denver 
2480  West  70th  Avenue 
Denver  21,  Colorado 

Mr.  James  W.  Sslassi,  President 
Turbocraft,  Inc. 

1946  S.  Myrtle  Avenue 
Monrovia,  California 

Universal  Match  Corporation 
Avionics  Department  Technical  Library 
4407  Cook  Avenue 
St.  Louis  13,  Missouri 

Waukesha  Bearings  Corporation 

p.  0.  Box  346 

Waukesha,  Wisconsin 

Attn:  Mr.  J.  M.  Gruber,  Ch.  Engr. 

Mr.  John  Boyd 

Westinghouse  Electric  Corporation 

Raaearch  Laboratories 

East  Pittsburg,  Pennsylvania 

Mr.  H.  Walter 
Director  of  Research 
Worthington  Corporation 
Harrison,  New  Jersey 

Dr.  W.  A.  Gross 
Ampex  Corporation 
934  Charter  Blvd. 

Redwood  City,  California 

profassor  J.  Modray 

Deparmaot  of  Machanlcal  Englnaarlnt 

Union  Cellaga 

gcbanactady  8,  t®*** 
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